Superconducting gap variations induced by structural supermodulation in BSCCO 
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We discuss the possibility that the strain field introduced by the structural supermodulation in 
Bi-2212 and certain other cuprate materials may modulate the superconducting pairing interaction. 
We calculate the amplitude of this effect, visible in scanning tunneling spectroscopy experiments, 
and thereby relate a change in the local superconducting gap with the change in the local dopant 
displacements induced by the supermodulation. In principle, since this modulation is periodic, 
sufficiently accurate x-ray measurements or ab initio calculations should enable one to determine 
which atomic displacements enhance pairing and therefore T c . 
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In the twenty years since the discovery of the high-T c 
cuprates, an intense theoretical and experimental effort 
has not reached a consensus on the origin of the phe- 
nomenon, although much has been learned. Progress 
has been slow due both to the complexity of the ma- 
terials, and the difficulty of the theoretical problem of 
strong electronic correlations. One advantage researchers 
in the high-T c field have is a set of new experimental 
tools which provide local information on the electronic 
state. Chief among these new methods is scanning tun- 
neling spectroscopy (STS), which is providing fascinat- 
ing new insights and forcing new ways of thinking about 
the high-T c _problem by yielding an unprecedented level 
of detail [H, H, & Hi HI- Traditionally, such atomic scale 
information has been considered irrelevant to the phe- 
nomenon of superconductivity, since rapid oscillations of 
pair wavefunctions on length scales smaller than the co- 
herence length £ are integrated out in the conventional 
pairing theory, where it is assumed that £ ^> a, with a the 
lattice spacing. In the cuprate superconductors, on the 
other hand, £ is much smaller than in conventional ma- 
terials, approaching a itself, and it is conceivable that a 
new type of approach accounting for atomic scale physics 
will be required to solve this problem. 

This possibility was highlighted recently when McElroy 
et aJ. [6] discovered the positive correlation of the posi- 
tions of dopant atoms in superconducting Bi-2212 with 
the local energy gap. It was argued that these atoms were 
in fact interstitial O atoms @, probably located between 
the BiO and SrO layers [7]. Subsequently, Nunner et aZ.Q 
showed that a good fit to several independent measured 
STS correlations could be understood if one assumed that 
the dopants, in addition to delivering holes to the Cu02 
plane and inducing a screened Coulomb potential, also 
modulate the local pair interactions. The idea is that 
each dopant distorts the lattice around it in such a way 



as to modify the effective electronic structure, character- 
ized by hopping integrals t, t', exchange constants J, 
electron-phonon couplings A, resulting in a modified local 
pairing interaction between electrons as welljfl. flOl . Phe- 
nomenological fits to the data of McElroy et ai.p then 
led to the conclusion that a substantial modulation of the 
pair interaction on an atomic scale is present in the ef- 
fective low- energy Hamiltonian of the disordered BSCCO 
material HEi 

More recently, Slezak et al.[l4| performed an exten- 
sive STS study of the supermodulation (SM) in BSCCO. 
From these local measurements it was found that the gap 
is modulated with the SM phase (j) SM with an amplitude 
of order 10% of its average value in near-optimally doped 
samples, providing a remarkable direct quantitative link 
between atomic displacements in the unit cell and the 
local superconducting gap A. In other words, the SM 
induces a pair-density wave in Bi-2212. 

In this work we extend the model of Ref. 8 to in- 
clude effects of the SM strain field. In a similar spirit, 
we assume that the local atomic displacements caused 
by the SM generate an additional, sinusoidal modula- 
tion of g, which we refer to as a g-wave. The su- 
permodulation presumably originates from a mismatch 
between the insulating layers and the preferred bond 
lengths of perovskite crystal. It can be characterized by 
a wavevector qsM, has a wavelength of approximately 
^sm — Zn/qsM ~ 26 A « 4.8 unit cells, runs along the 
a-axis, and leads to displacements of the atomic posi- 
tions of up to 0.4 A. It is believed that one of the main 
effects of the SM is to modulate the distance between the 
apical oxygen and the Cu0 2 plane 0, [H, [HI- 

For sim- 
plicity, we neglect the periodic changes in other terms in 
the Hamiltonian, such as the electron hopping. We will 
show that the amplitude of this g-wave can in principle 
be determined by comparison with experiment, given the 
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assumptions already specified. In this case, one should 
be able to relate the change in the pair potential Sg to 
the atomic displacements caused by the SM, informa- 
tion which should be available with sufficiently precise 
x-ray data, or from ab initio calculations. The philoso- 
phy here is similar to that of Nunner et al.[8], but the 
measurement of changes in the local gap caused by the 
supermodulation has the advantage that the associated 
periodic displacements should be easier to determine em- 
pirically than in the case of the random O dopants. As 
discussed below, current x-ray data on Bi-2212 are not 
yet able to answer this question, but there appears to be 
no fundamental obstacle to improving the precision so as 
to be able to achieve this goal; they may then be able to 
remove any remaining ambiguity as to the microscopic 
origin of the modulated pairing in this material. 

In the following, we use the conventional d-wave BCS 
Hamiltonian defined on a 2D square lattice 

h= J2 *«4cja+^K-/*)4cfa+5Z ( A «44 +Hx ) ' 

(ij)" ia (ij) 

(1) 

where in the first term we include nearest t and next- 
nearest t' = — 0.3t neighbor hopping. We set the chemi- 
cal potential /i = — t in order to model the Fermi surface 
of BSCCO near optimal doping, and denotes sum- 

mation over neighboring lattice sites i and j. Disorder of 
the usual screened Coulomb type is included in the impu- 
rity potential V { = V fi where = ^ s exp(-r is /A)/r is , 
where r is is the distance from a defect s to the lattice site 
i in the plane. Distance (energy) is measured in units of 
\/2a (£), where a is the Cu-Cu distance, and the calcula- 
tions are done at T = 0. Note that the particular form 
of fi is merely a convenient way to vary the range of the 
impurity potential landscape through the parameter A. 
The d-wave order parameter A^ = gij(c^Cji — Cj|C^)/2 
is determined self-consistently via 

A ij = ~Y ( u n(i)vn(j) + v n (i)u n (j)) tanh(^), (2) 

n 

Here, {E n , u n , v n } is the eigensystem resulting from di- 
agonalization of the Bogoliubov-de Gennes (BdG) equa- 
tions associated with Eq.(pQ). The pairing interaction g^ 
is assumed to vary in space relative to its homogeneous 
value go as 

9ij =9o + Sg imp (fi + fj)/2 + SgsM cos(0f M ), (3) 

where Sgi mp and SgsM are the amplitudes of the dopant 
and g-wave modulation, respectively, and z, j are nearest 
neighbors. In the homogeneous case g^ = go we choose 
go = 1.16t giving A^ = O.lt. In the inhomogeneous 
case we make sure to adjust go, m order to maintain the 
same average gap as in the corresponding homogeneous 
system. The third (g-wave) term in Eq. (|3]) is taken to 
vary sinusoidally with the SM phase (/>f M , a phase vari- 
able associated with the structural supermodulation as 



determined in Ref. [3- In the following, we use the ex- 
perimentally determined <pf M as input in Eq. (j3]), and 
ignore for simplicity the conventional potential (r^ chan- 
nel) associated with each dopant, and consider disorder 
only in the pair (n) channel. The goal will be to ad- 
just the amplitude SgsM to obtain better quantitative 
agreement with the experiments of Ref. LL4|. 

That a nonzero value of 5g$M is required is illus- 
trated in Fig. [T] which shows results obtained with 
SgsM = 0[12]. The input to the theory is, in a field 
of view (FOV) of 49nm x 49nm, the measured conduc- 
tance at -960 meV (a), known to track the positions of the 
dopant oxygens. The corresponding gapmap obtained by 
McElroy et al.[6] is shown in (b). In Fig. [He) we dis- 
play the impurity potential generated by assuming that 
each of these dopants provides a potential centered on the 
bright spot positions of (a) which decays as a screened 
exponential in the pairing channel, as described above. 
The similarity to (a) is manifest. The experimental FOV 
is modeled as a 90 x 90 lattice system rotated 45 degrees 




FIG. 1: (Color online) (a) Experimental dl/dV map [arb. 
units] at -960 meV of an optimally doped BSCCO sample^], 
(b) Experimental gapmap [meV] in the same region as (a) 
at T = 4K. (c) The theoretical impurity potential extracted 
from (a) assuming a distance from the Cu02 plane r z = 0.5 
and A = 0.5. (d) The gapmap resulting from using (c) as the 
pairing potential in the BdG equations with Sgi mp = 1.5t. (e) 
and (f) display the ^-histogram and A-histogram versus (f) SM , 
respectively. The model results (c-f) are for SgsM = 0.0. 



with respect to the 3.83 A Cu-Cu bonds. Therefore our 
system consists of 2 x 90 x 90 sites and is aligned with 
the experimental FOV. The gapmap computed from the 
coherence peak-to-peak distance in the theoretical local 
density of states (LDOS) using (c) as input pair potential 
8gimp is shown Fig. QJd). Gapmaps reasonably consis- 
tent with experiment are found for 5gi mp ~ 1.5£. 

Although various correlations among O positions and 
the LDOS p{E) are successfully reproduced by the Nun- 
ner et al [8] approach, one deficiency appears upon 
closer examination of the gapmaps Fig. []Jb) and (d). 
The experimental result (Fig. [U(b)) contain nearly verti- 
cal linear striated modulations visible to the eye, which 
match the corresponding STM topograph and are there- 
fore caused by the SM. These are not obviously manifest 
in the theoretical gapmap presented in Fig. [T] (d). In 
fact, weak correlations of the dopant positions with the 
supermodulation phase are indeed present in the data, as 
seen in Fig. [U(e-f) and lead, for the disorder parameters 
which reproduced the gapmap statistics, to a small net 
modulation of the gap with <p SM . 

In Fig. HJe), it is seen that the experimental dopant 
distribution is correlated with <p SM at and 180°, lead- 
ing to a peaks in the coupling g vs <p SM at these phases 
in our model. Fig. [Iff) shows a histogram of the gap 
in Fig. QJd) versus SM phase (j) SM . Because of the ef- 
fective smearing of the order parameter response over 
scales of order the coherence length £ ~ Asm, the 180° 
peak is partially wiped out and a resulting weak 360°- 
periodic modulation of 4-5% of A versus (j) SM remains, 
qualitatively similar to the experiments [141]. Increasing 




FIG. 2: (Color online) For clarity we show here the self- 
consistent results for the superconducting gap A resulting 
from a pure g-wave (Sgimp = and SgsM = 0.08£): (a) 
real-space gapmap, (b) conventional gap histogram, (c) g- 
histogram versus <j) SM , and (d) A-histogram versus (f) SM . 
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FIG. 3: (Color online) (a,b) Real-space gapmaps for Sgimp = 
1.5t, and SgsM — 0.06£ (a) and SgsM = 0.14£ (b). For the 
parameters in (a) we show is (c) the ^-histogram versus <j) SM , 
and in (d) the A-histogram versus <j) SM . 



the phenomenological amplitude of the dopant potential 
leads to overall fluctuations of the gap amplitude which 
are much too large compared to experiment. Therefore, 
in order to generate a 10% variation of A versus (j) SM , 
and maintain quantitative agreement with the spatial ex- 
tent and overall amplitude of the gap variations evident 
from the experimental gapmaps, we need to include a 
nonzero 5gsM- This is fully consistent with the notion^ 
that the reason the oxygens modulate the pairing is due 
to local distortions of the lattice, and the SM should pro- 
duce similar effects on its own. Hence, in the following 
we investigate the possibility that a g-wave is present in 
the system in addition to the dopant t\ disorder. 

In Fig. [21 we first show the effects of a pure g-wave 
without r\ disorder, i.e. Sgi mp = 0. From Fig. EJa), 
which displays the gap in real-space, it is seen that the 
SM agrees well to that observed in the corresponding ex- 
perimental FOV by Slezak et al. [14| as it should per con- 
struction. The gap histogram is shown in Fig. [21(b), and 
in Fig. [2](c) (Fig. [21(d)) we show the histograms for the 
input g-wave (self-consistent A) versus SM phase <p SM , 
both of which trivially exhibit a sinusoidal dependence. 
Fig. [3] displays typical results when including both Sgi mp 
and SgsM into the simulation. Here, we used the same 
disorder parameters as in Andersen et al.[12] to obtain 
the gapmaps (Fig. [Hc,d)), but with SgsM — 0.06£ (a) 
and SgsM — 0.14£ (b). In both these plots the striated 
gap modulations are evident. For the parameters used 
in Fig. [3](a), we show in EJc) and [31(d) the g-histogram 
and A-histogram, respectively. As seen, the oscillation 
of g versus (j) SM of roughly sinusoidal form is still found, 
with amplitude close to the input SgsM- The gap re- 
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sponse has a corresponding amplitude of about 10% of 
its average value in agreement with the measurements by 
Slezak et al. [13] • Therefore, within the modulated pair- 
ing scenario, an experimental gap modulation amplitude 
of approximately 10%, is caused by a SM-induced g-wave 
of similar size relative to the background pairing strength 
go. Note that after including the SM phase, it is Fig. [3ja) 
that should be compared to the experimental gapmap in 
Fig. [Kb). 

We would now like to extract which type of atomic 
displacements are present in the SM and associated with 
the enhanced pairing at and 360 degrees SM phase. In 
principle, this information should be available from care- 
ful x-ray diffraction data, but this is complicated by the 
existence of deviations of this system from stoichiometry, 
both through Bi/Sr substitutions and oxygen intersti- 
tials, which together determine the incommensurabi lity 
of the true system. Early x-ray analysis TH, 3, 2(J 21] 
suggested a weak correlation of the interstitial oxygen 
position with the SM, but as remarked above this cor- 
relation does not appear to be the most significant one 
present in the STM gapmaps. 

More recently, progress in the analysis of incommensu- 
rate systems 
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231 ] has been made, and new short-range 
ordered structures have been identified [24], but the var- 
ious studies disagree at essential points[14], and it does 
not appear to be possible with present data to identify 
the actual displacements of atoms near the SM maxi- 
mum or minimum with sufficient precision to eventually 
draw conclusions regarding the microscopic origin of the 
pairing modulations. There does appear to be a statis- 
tically significant correlation of the z coordinate of the 
apical oxygen relative to the Cu02 plane with the SM 
phase [25], but it anti-correlates with the local gap [26]. 
This behavior is the opposite of that which might be ex- 
pected on the basis of the analysis by Pavarini et al. [13] , 
who pointed out an empirical correlation between apical 
O coordinate z and T c . This may point to the primacy 
of other atomic displacements, or suggest that T c itself 
is less related to in-plane pairing strength and more to 
inter layer couplings. In any case, our analysis should pro- 
vide incentive for a repeated assault on the x-ray analysis 
of this compound, with the hope of eventually providing 
direct local information on the origin of the pairing. 
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